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Abstract 

The low-energy effective theories for the two- and three-flavor color- 
superconductors arising in the high density limit of QCD are discussed. Using 
an effective field theory to describe quarks near the fermi surface, we compute 
the masses of the pseudo-Goldstone bosons that dominate the low-momentum 
dynamics of these systems. 
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I. INTRODUCTION 



Recent developments have reinvigorated efforts to understand QCD at very high baryon 
density [|1^. For special combinations of quark colors and flavors it is likely that a 
superconducting gap breaks color and flavor symmetries in interesting ways. Although the 
symmetry breaking is nonperturbative, it occurs when QCD is weakly coupled, and therefore 
perturbative QCD (pQCD) can be used to derive properties of the superconducting phase. 
One can conceive of scenarios in nature where it may be important to understand the 
behavior of QCD at high density: for instance, in neutron stars and less likely, in high 
energy heavy ion collisions. 

Below the weak scale the standard model has the exact local gauge symmetries S'f/(3)c® 
6^^(1)6111 which describe the strong and electromagnetic interactions. In addition there is the 
exact global symmetry U{1)b corresponding to the conservation of baryon number and the 
approximate global symmetries SU{Nf)L (g) SU{Nf)ji U{1)a for a theory with Nf = 2,3 
light quarks. These global symmetries are broken by the quark masses and in addition 
the U{1)a symmetry is also broken by the strong anomaly. At extremely high density, the 
contribution from the anomaly is suppressed by powers of the chemical potential, /i, and 
U{1)a is broken only by the light quark masses. When a color-superconductor forms via 
the spontaneous breaking of color and flavor symmetries, there will be pseudo-Goldstone 
bosons that contribute to or determine the very low energy dynamics of such a system. For 
sufficiently low energies it is clear that an effective field theory description of these dynamics 
can be constructed, and will prove useful in computing contributions to observables from 
the far-infrared region of the theory. In recent work by Son and Stephanov |]l3l and by 
Casalbuoni and Gatto Ijl^ the masses and decay constants of the pseudo-Goldstone bosons 
in the Nf = 3 color-flavor scenario were computed in the large-/i limit. It was found that 
the masses become independent of fi, while the decay constants depend linearly on fi. Later 
work in [T^ and |16| claim that the masses actually vanish in the large-/! limit and are 
proportional to ~ A^//i^ log(A//i). We agree with this later claim |1^,|T5|JT^ and, through 
the use of a hierarchy of effective theories, compute the leading contribution to the meson 
masses and decay constants. 



II. THREE FLAVORS 

If we assume that the masses of the up, down and strange quarks are much smaller than 
the scale associated with the formation and dynamics of the color-superconducting state (the 
gap A), then it is appropriate to consider a theory with three flavors of massless quarks, 
and include mass effects in perturbation theory. In the limit of high densities the attraction 
leading to the gap is given by one-gluon exchange that is attractive in the color 3 channel. 
It has been argued that the most favorable state for the system is one in which there is a 
formation of the "color-flavor locked" condensate 

(Ka^^j) = = A6"^^6,,,6,, (l) 

{a, (3,.. are color indices, a, 6,. .are flavor and are spin indices) resulting in the 

symmetry breaking pattern S't/(3)c ® SU{3)l ® SU{3)r ® U{1)a ® U{1)b ® ?7(l)cm ^ 
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SU{3)c+L+R ® f^(l)cm- The value of the gap A was computed in PJTO|J^,pl]] in the high 
density hmit and does not have the usual BCS form but is instead given by 

where c is a constant of order unity not yet computed. For lower densities, where pQCD 
does not apply, the same symmetry breaking pattern was shown to occur by assuming that 
all quark interactions are effectively short ranged. Throughout this work we assume that A 
is a constant, independent of energy and momentum. 

As we are interested in modes close to the Fermi surface where |p| ~ ;U ^ A we can 
start by considering the dynamics of quarks and gluons in the ungapped system, where the 
dynamics are determined by the 3 + 1 dimensional action, 



S3+1 = I dt d X 



(3) 



It is convenient to project onto the positive and negative energy states, ip^ and ijj- respec- 
tively, and then eliminate ip^ using the equations of motion 0, with 

iIj = ^+ + ij. , ^± = , P± = ^ (1 ± ToT^n'^) , (4) 

where n = p/|p|. This procedure is similar to that used in the construction of heavy quark 
effective theory (HQET) |jl9|. Writing S^+i in terms of the mode expansion for ipj^ and 
working in spherical coordinates, the action describing the dynamics of the modes near the 
fermi surface, iSs+i, is 

'^3+1 = — , (5) 

TT 

where Si+i is the action of a 1 + 1 dimensional field theory. The two-component quark field 
= Xa(-^!^)J^) depends on an energy <C /i, a momentum k <^ fi, and a unit vector 
pointing toward the fermi surface n. As the anti-quarks have been integrated out of the 
theory at the scale /i corresponding to the top of the fermi surface, the effective field theory 
described by Si+i will be an expansion in terms of E/fi and k/fi, as outlined in 0. 

Perturbative computations around the superconducting state can be performed by adding 
and subtracting a quark gap term in the QCD lagrangian. The condition that the subtracted 
gap term does not contribute at each order in perturbation theory is equivalent to the gap 
equation. The gaps for the positive and negative energy states will in general be different, 
and in terms of the ip± fields we have an additional contribution to the lagrange density of 
the form 

= |e.^,6'^^Vr+C< + f e^^/e'^^'V^r.-CV^f- + h.c. . (6) 

The anti-gap. A, has not been computed at this point in time, but a recent discussion can 
be found in . 

A further simplification can be made by writing the quark fields, iii terms of the mass 
eigenstates of the condensate (neglecting the quark masses) 
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Xa 

A=l 



in 



(7) 



where with A= 1, .., 8 are the Gell-Mann matrices and = y 2/3 Is. After ehminating 
ip- by the equations of motion, the part of the leading order action that does not depend 
upon the quark masses is 



c(0) 
'-'l+l 



^ rdndEdk 

A=l 



A-n (27r) 



A 



(8) 



where A"^ = A for A = 1, ..8, while A^ = — 2A. Interactions between the ^ind the gauge 
fields have not been shown. The ellipses denote operators that are suppressed by powers of 
/i. As we are assuming that the quark masses are small compared to the scales associated 
with formation of the superconducting state, we can treat the quark masses in perturbation 
theory. The leading order contributions from the quark masses are described by the action 



<.(m2) 
'-'l+l 



E 



dn dE dk 



Tr [A^A^tA^A^ 



A 



16/i2 



C^xI-„ + xflC^xf-n) 3^^^ + h.c. 



where 



y 



AB 



Tr 



Tr 



Tr 



(9) 



(10) 



and where the L, i? subscripts on x"^ denotes the helicity/chirality state. The appearance 
of explicit factors of 1//^ associated with the mass terms is no surprise, in fact, the first 
term in Eq. follows naturally from the expansion of £" = ^J'p^ + = p + /{2p) + .. 
for p fi ^ m. Contributions from more insertions of the light quark mass matrix or 
from higher derivative operators are suppressed by powers of fi. The actions in Eq. (^) and 
Eq. (1^) describe the dynamics of modes near the fermi surface. Contributions to observables 
arising from modes far from the fermi surface are suppressed by powers of /x, and enter 
through the higher dimension operators that we have not shown. An important point is 
that this lagrange density involves only analytic functions of the light quark masses and the 
gap A. For very low-energy dynamics of the system |p| ^ A it is appropriate to construct 
an effective field theory for the pseudo-Goldstone modes alone. This lagrange density will 
involve analytic functions of the light quark masses, but it will have non-analytic dependence 
on A. This nonanalytic dependence can be computed from the effective theory describing 
the momentum region A -C |p| <^ /i, as described above. This is in direct analogy with 
the nonanalytic contributions to observables in the light meson sector, such as terms of the 
form log (^m^/A^^ or ^m^. 

For momenta much below the gap, p -C A, the relevant degrees of freedom are the nine 
pseudo-Goldstone bosons resulting from symmetry breaking due to the condensate. They 
are described by the fields 



V 



(11) 



4 



with 







71+ 


K+ 


M = 










I K- 







(12) 



These fields and the quark mass matrix transform under SU (3) l x SU (3)^? x U{1) b x U{1)a 
as 



(13) 



where (3 is the U{1)a phase. The interactions of these mesonic excitations will have an 
expansion around zero momentum described by 



Tr 



V8 



f2 



f2 

+ 4f doVdoV* 



f2 
|2il 



+ Ai (Tt[M^ ^]Tt[M^ ^]V* + h.c.) + A2 (Tr [M^ ^M^ ^]V* + h.c. 
+ A3 Tr [ TW^ S ] Tr[ M St ] + Tr [MM^ ] + ■■■ , 



(14) 



where the ellipses denote operators suppressed by powers of p/A and Ai/A (for a recent 
discussion of this see |1^). The velocities of the octet and singlet bosons are vg and vi 
respectively. The coefficients appearing in Eq.(|l^ are determined by matching S-matrix 
elements in the theory above A with those in the theory below A, determined from Eqs.(|^), 
(^) and Eq. ( p^lj) respectively. The form of the mass terms is exactly the same as the second 



order terms appearing in the Lagrange density written down by Gasser and Leutwyler [20 



A. Decay Constants 



The decay constants fs and /i can be found by computing the Debye and Meissner 
masses of fictitious gauge bosons coupled to currents in both pQCD ||13| and the low-energy 
effective field theory, Eq. (JjL^j. In the low-energy regime the two calculations must produce 
identical results, thereby determining /§ and f\ order by order in the chiral, l//i and 
expansions. 

Gauging the U{1)a axial current leads to = + ieW^Qi, where Qi is the axial 
charge operator, Qi^l = +1^l, Qi'^r = -l^ij, QiV = -W and Qit = 0. It is then 
straightforward to show that the mass of the fields at leading order in the low-energy 
effective field theory are 



A£ = 2 /2 



(15) 



Computing these masses in pQCD through Eq. (||) one finds that the Debye mass is given 
by the two graphs shown in Fig. (|I]) with the result 



AC 



QCD 



and is essentially identical to the standard many-body calculations described in |2T 
order to reproduce this with the effective theory 



(16) 
. In 
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(a) (b) 

FIG. 1. One loop diagrams that give the leading contribution to the meson decay constants. 



fi 



4 7]-2 



(17) 



which agrees with To compute the speed of the Goldstone mode, we need to compute 
the Meissner mass in both theories. In addition to the diagrams of Fig. (|l]) (with different 
couplings) there is a contribution from a counterterm. The sum of diagrams yields a vanish- 
ing Meissner mass in the normal phase and a non-zero mass in the superconducting phase, 
see The speed is found to be |vip = i, in agreement with |T^. As discussed in the 
analogous calculation for the baryon number Goldstone boson is identical at leading order, 

12 _ U, 12 _ 1 
|V//| - |Vl| - 3. 

To determine the decay constant of the octet pseudo-Goldstone bosons, fs, the Debye 
and Meissner masses of fictitious gauge bosons associated with octet currents are computed. 
We find 



fs 



/i^ 21-81og[2] 
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(18) 



which differs by a factor of 2 from Ref. [|T^ after differences due to the definition of fs are 
taken into account. The pion velocity is found to be Ivgp = | which agrees with Ref. |]13 . 
Therefore, all Goldstone modes have the same speed at leading order in the expansion. 

As the decay constants scale like / ~ /i, it is apparent that the dynamics of the pseudo- 
Goldstone bosons do not receive significant contributions from loops in the low-energy ef- 
fective theory of Eq. ([T^). The naive size of the counterterms is set by 1/A, while the 
contribution from loops is set by Therefore, once the coefficients in the low energy 

effective theory have been determined, only tree-level diagrams need to be considered. 
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(a) (b) (c) 

FIG. 2. One loop diagrams that give the leading contribution to the meson masses. Graph (a) 
contributes only to A4 which does not impact the dynamics of the pseudo-Goldstone modes, while 
graphs (b) and (c) contribute to both Ai and A2. The crossed circle denotes an insertion of the 
quark mass operators in the effective field theory defined in Eq. Ij^. 



B. Meson Masses 



The coefficients in Eq. (|14D can be determined by matching the change in the ground 
state energy due to the quark masses in the high (Eqs. (H) and (^) and low energy theories 
(Eq. (0)). In the low energy theory, the change in energy density can be easily determined 
from Eq. (|T^ by setting S = VI^ = /a, to yield 



5£ = ((Tr [M]f + h.c.) + A2 (Tr 



+ h.c. 



+ AgTr [M] Tr \M^ 



(19) 



The operator with coefficient ^4 does not contribute to the dynamics of the pseudo- 
Goldstone modes, and we do not calculate it. 

Computation of the energy density in the 3 + 1 dimensional high energy theory can easily 
be done, by noting that the action of the 3 + 1 dimensional theory is a factor of /x^/tt times 
the action of the 1 + 1 dimensional theory. Thus the energy density in the 3 + 1 dimensional 
theory is /i^/vr times the energy density computed in the 1 + 1 dimensional theory. We use 
dimensional regularization and minimal subtraction to define divergent integrals that occur 
at loop level in the 1 + 1 dimensional theory. The shift in the vacuum energy due to the 
light quark masses results from the tadpole diagrams shown in Fig. (Q), where the vertex 
arises from Eq. (Rl). At leading order we find 
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^log(^) ( (Tr[^])^ - Tt\m'] ) + h.c. 



47r2 ^'A' 



(20) 



where A is the renormahzation scale and we have only shown the term nonanalytic in A/ A. 



The form of our expression agrees with the results of [|I2],|I5|,|I6| . The explicit dependence on 
A shown in Eq. (|20|) is absorbed by an equal but opposite A dependent counterterm that, for 
A ~ yU generates a shift in energy ~ AA/An"^. This counterterm contribution is suppressed 
compared to the contribution in Eq.(^) by the large log (A/yu) factor and will be neglected. 
Matching this result with the corresponding shift in energy computed in the effective theory 
we find 



-A, 



^log — 



A , A, 







(21) 



The meson masses at leading order are found by expanding Eq.(|l4D to second order in 
the meson fields. The charged meson masses are 

= ^ (rriu + md)ms , m\+ = ^ {rriu + ms)md 
Is Is 



m 



8A 

li 



{rrid + ms)mu 



221 



and the neutral meson mass matrix is 

2 _ 8^ 
33-72" 

IQA 



m 



■88 



8A_ 



[ msiniu + rrid) + ^rriumd ] 



m 



11 



IQA 



msiniu + TTid) + rriumd ] , m 



8AV2 



13 



/i/s 



rrid rris 



m 



18 



msiniu + md) - 2mumd ] , m] 



38 



8A 



nid ms 



(23) 



V^/l/8 

In the limit of zero mixing mn = m^/, 77133 = i^-k" and mgg = nir^. The eigenvectors of this 
matrix are close to the quark flavor eigenstates for values of the quark masses consistent 
with standard chiral perturbation theory ||20|| . 

It is interesting to consider mass relations between the mesons. At this order in pertur- 
bation theory there is a relation between the meson masses without any assumption about 
the light quark mass hierarchy, 



2 I 2 I 2 
m^, + m^o + 



2 ^ 
^ f? 



m^ 



+ ni\+ + m\o 



(24) 



where m^, m^/ and m^o are the mass eigenvalues found by diagonalizing the neutral meson 
mass matrix. Further, if one assumes mu^d/^s ^ 1 and neglects such terms. 



m 



nid 



m 



m\o 
2 

m + 



(25) 



and hence + m|-o 



md + mu m;+ + m„ 

m^+. It is also clear that there is an inverse mass hierarchy, e.g. 
mT,+ > mK+ > m^o. 

In addition to corrections to the coefficients Ai and A2 arising at higher orders in l/n 
and as there will be contributions to the ^3, such as those shown in Fig. (|^). The operator 
in the effective theory appearing in Fig. (|^) results from a two loop diagram in QCD where 
two of the propagators are far off-shell, resulting in an effectively local vertex for momenta 
much less than /i. 
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(a) 



(b) 



FIG. 3. Some of the two loop diagrams that contribute to A3. The crossed circle denotes an 
insertion of a quark mass operator involving the gluon field that arises at higher orders in the 
effective field theory expansion. The curly lines denote gluons. 



III. TWO FLAVORS 

In the two flavor case (three colors) the most favored condensate is 

{na^<b,) = ina^Kt,) = Ae^^^e^.e,, (26) 

which breaks SU{3)c SU{2)l O SU{2)r U{1)a ® U{1)b O f/(l)em down to SU{2)^ O 
SU{2)l ® SU{2)ii ® f/(l)em ® U{1)^, where f/(l)em and f/(l) are linear combinations of 
electromagnetism, baryon number and the eighth gluon, PQ- At energies below the gap 
the dynamical degrees of freedom are three massless gluons, two ungapped quarks and the 
pseudo-Goldstone boson associated with spontaneous breaking of U{1)a- In contrast to the 
three flavor case, there is no pseudo-Goldstone boson associated with the breaking of baryon 
number. The Goldstone field describing the 77' meson related to the breaking of U{1)a can 
be parametrized as 

y = g^2r,'//v ^ v^e-'^^V . (27) 

The effective lagrangian describing the dynamics of the rj' at leading order in the d/A and 
m/A expansion is 
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C^-^( DqVDoV* - \v\'^DiVDiV* ) + B [det (M) V + h.c. ] . (28) 
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Using the diagonal basis for the gapped quark fields with a, i = 1, 2 

*" = ;^EM>^ (29) 

where are PauU matrices for A = 1, 2, 3 and ~ I2, the gaps for the four gapped quarks 
are — A4 = = A for A = 1,2,3, and the others remain ungapped. Calculation of the 
decay constant and vacuum energy shift due to the light quark masses is similar to that 
performed in the three flavor case and gives 

/^.^ , M-i . B = -^M^) (30) 

and therefore an rj' mass of 

2 SBmumd 4AA A 
"^r,' = 72 = ^log(— )m„md , (31) 

Jr]' A* 

where we have neglected the contribution from the local counterterm to the r]' mass. 

The ungapped quarks interact with themselves at leading order via the exchange of a 
massive gluon, inducing a four-quark operator for scales below gsfJ.- The coefficient of this 
operator is independent of the strong coupling constant ag (which arises from a cancellation 
between the couplings and the gluon mass) and scales like 1///^. This interaction is repulsive 
since the two ungapped quarks are in a 6 of color and further the interaction vanishes in the 
high density limit. Thus we do not expect condensation of the "green" quarks in the high 
density limit. It is interesting to note that the rj' does not couple directly to the ungapped 
quarks since the rj' field is an excitation of a condensate involving the "red" and "blue" 
colors only. Thus the axial coupling constant describing this interaction is suppressed by 
powers of ag and consequently the tensor force between the ungapped quarks in the low 
energy theory arising from the exchange of a single rj' is suppressed by a^. 

Another interesting aspect of the two-flavor case is the presence of an unbroken pure 
SU{2)c gauge theory. The gluons associated with this gauge group do not interact with the 
ungapped quarks. The conflnement scale of this theory can be estimated by assuming the 
only modification to the evolution of the strong coupling arises from the particle content. 
This provides an estimate of the scale at which the theory becomes strongly coupled, 

__7_ 

;^conf. _ ^conf. ( ) e^"(^-^) , (32) 

V^QCD/ 

where fJ^f^^' are these scales in the two and three color theories respectively. For large A, 
this scale is much lower than Aqcd, but for reasonable values of A, /i2°°^' ~ /x™'^^'. This 
suggests that pure Yang-Mills glueballs will appear in the low energy theory with masses of 
order ~ Aqcd- 
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IV. CONCLUSION 



We have examined the high density hmit of QCD where there are two and three flavors 
of "hght" quarks below the scale relevant to the formation of a color superconducting state. 
Using an effective field theory to describe quark modes near the Fermi surface we have 
determined the decay constants and masses of the pseudo-Goldstone bosons that arise in 
each theory at leading order in the expansion. The masses of these pseudo-Goldstone 
modes vanish in the high density limit. In order to determine the behavior of these systems 
at a moderate density the subleading corrections (e.g. ag) need to be determined, 

including the contributions from instantons ITH . 
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